arXiv:1502.05878v2 [gr-qc] 6 May 2015 


Prepared for submission to JCAP 


Attracted to de Sitter: 
cosmology of the linear Horndeski 
models 


Prado Martm-Moruno, Nelson J. Nunes and Francisco 
S. N. Lobo 

Institute de Astrofisica e Ciencias do Espago, Faculdade de Ciencias da Universidade 
de Lisboa, Campo Grande, PT1749-016 Lisboa, Portugal 

E-mail: pmmoruno@fc.ul.pt, njnunes@fc.ul.pt, fslobo@fc.ul.pt 

Abstract. We consider Horndeski cosmological models, with a minisuperspace La- 
grangian linear in the held derivative, that are able to screen any vacuum energy and 
material content leading to a spatially hat de Sitter vacuum hxed by the theory itself. 
Furthermore, we investigate particular models with a cosmic evolution independent of 
the material content and use them to understand the general characteristics of this 
framework. We also consider more realistic models, which we denote the “term-by- 
term” and “tripod” models, focusing attention on cases in which the critical point is 
indeed an attractor solution and the cosmological history is of particular interest. 
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1 Introduction 

Scalar-tensor theories of gravity have recently received a renewed interest in relation 
with the problem of describing the current accelerated expansion of our Universe. In¬ 
deed, scalar helds are ubiquitous in theories of high energy physics beyond the standard 
model and, in particular, are present in theories which include extra spatial dimensions, 
such as those derived from string theories. The initial simplest model, denoted by the 
quintessence held [1], consisted of a minimally coupled single scalar held self-interacting 
through a scalar potential, and with canonical kinetic terms. An interesting model, de¬ 
scribed as coupled quintessence, which analysed the consequences of coupling the held 
responsible for acceleration to matter helds [2] , was rapidly extended to consider non- 
canonical kinetic terms, with higher powers of the held’s velocity, denoted fc—essence 
[3], and non-minimal couplings to gravity. Scalar helds with Lagrangians presenting 
nonlinear kinetic interactions and possessing shift symmetries, known as galileons, were 
found to appear in the decoupling limit of the DGP models [4]. However, the “self- 
accelerating” solution in the DGP model is plagued by a ghost instability, rendering 
the solution untenable. This fact motivated the proposal of an infrared modihcation 
of General Relativity (GR) [5], involving an internal “galilean” invariance, where the 
symmetry constrained the structure of the ehective Lagrangian without introducing 
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ghosts. These galileon helds have second order equations of motion and lead to the 
Vainshtein screening mechanism at astrophysical scales [6]. Galileon scalar helds also 
appear in the decoupling limit of the ghost-free theory of massive gravity [7, 8]. The 
understanding of this theory [9, 10] as well as of the intimately related [11] bigravity 
theory [12] , free from the Boulware-Desser ghost [13] has also recently been extensively 
explored (we refer the reader to Ref. [14] for a recent review). 

In fact, given the large number of models, the question arises how one should 
study and compare them in a unihed manner, and determine which if any is the origin 
of cosmic acceleration. A particularly useful tool in this direction is the realisation 
that the covariantized Galileons [15] (recovered in Ref. [16]) are special cases of the 
most general Lagrangian which leads to second order field equations [17], and there¬ 
fore yielding the model free of Ostrogradski instabilities. (This is, of course, also the 
case for the covariantization of the decoupling limit of massive gravity presented in 
Refs. [18, 19].) This Lagrangian was hrst written down by Horndeski in 1974 [20] and 
recently rediscovered by Deffayet et al. in 2011 [17] in an equivalent formulation [21]. 
An appealing theoretical property is the possibility that the held is able to self-tune 
considering some terms within the Horndeski Lagrangian (see Refs. [8, 22] for early 
works on degravitation). This can be used to partially alleviate the cosmological con¬ 
stant problem, which is not only based on the huge hne-tuning needed to cancel the 
value of the vacuum energy in particle physics to produce the value of the cosmological 
constant as inferred from cosmological observations, but also on the need to hne-tune 
again this difference after any phase transition. In particular these terms have been 
investigated with the aim of hnding a viable self-tuning mechanism, which screens 
the spacetime curvature from the net cosmological constant [14]. Horndeski models 
are able to self-tune any vacuum energy leading to a Minkowski vacuum space have, 
therefore, been obtained [23, 24]. 

Starting from the most general scalar-tensor theory with second order held equa¬ 
tions in four dimensions, the unique action that allows for the existence of a consistent 
self-tuning mechanism on FLRW backgrounds was established in Ref. [23]. It consists 
of a combination of only four base Lagrangians, denoted the Fab Four, with an intrigu¬ 
ing geometric structure dependent on the Ricci scalar, the Einstein tensor, the double 
dual of the Riemann tensor and the Gauss-Bonnet combination. It was shown that 
spacetime curvature can be screened from the net cosmological constant at any given 
moment, as the scalar held is allowed to break Poincare invariance on the self-tuning 
vacua, thereby evading the Weinberg no-go theorem [25]. This is only possible, how¬ 
ever, when there is a nonvanishing spatial curvature, a requirement that has raised 
instabilities related concerns in the literature [19] . The cosmology of the Fab-Four La¬ 
grangian, which also screens the contribution of any matter content, has been studied 
in Ref. [26]. By performing a phase plane analysis of the system, a number of hxed 
points for the system were obtained. In addition to obtaining inhationary solutions, 
conventional radiation/matter-like solutions were also found, but in regimes where the 
energy density is dominated by a cosmological constant, without the explicit forms of 
radiation or matter. An extended period of ‘matter domination’ was shown to exist, 
which opens up the possibility that cosmological structures can be generated, and thus 
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recovering a consistent cosmology before the Minkowski attractor is reached. While 
not yet a full theory of cosmological evolution, such a mechanism is intriguing, and 
leads to the question of whether realistic cosmologies that include this mechanism as 
well as the late time acceleration, can be constructed. 

Some of the terms appearing in the Fab Four models [23, 24] have been generalized 
[27]. More specihcally, it was shown that a nonlinear combination of purely kinetic 
gravity terms can give rise to an accelerating universe without the addition of extra 
propagating degrees of freedom on cosmological backgrounds, and exhibit self tuning to 
bring a large cosmological constant under control. However, concerns of the viability of 
the cosmologies described by this generalization have been raised in Ref. [28]. Indeed, 
these generalizations may be ghost free and stable, have attractor solutions in the 
past and future, and possess self tuning that solves the original cosmological constant 
problem, however, it was shown that they do not possess all of these qualities at the 
same time. It was also demonstrated that the screening is so powerful that it not only 
cancels the cosmological constant but also all other energy densities. Nevertheless, it 
must be noted that this is a general characteristic of self-tuning models, and it only 
entails a problem for models in which the self-tuning takes place so fast that it preclude 
the existence of a matter dominated phase before the accelerated expansion. 

More recently, Horndeski models which can screen any vacuum energy to a spatial 
flat de Sitter vacuum hxed by the theory itself have been investigated [29] . Thus, the 
most general minisuperspace Lagrangian for a scalar held leading to an equation of 
motion without higher than second derivatives and a spatially hat de Sitter critical 
point has been presented [29]. These models have the advantage of allowing to under¬ 
stand the current epoch of accelerating expansion as the result of the approach to a 
stable critical point while avoiding the cosmological constant problem. In this paper 
we investigate the cosmology of a family of these models which have a minisuperspace 
Lagrangian linear on the held derivatives. 

This article is outlined in the following manner: In Section 2, we present the 
model and include some expressions that will later be useful for our study. In Section 
3, we present some considerations regarding the dynamical problem and conclude that 
one needs to go beyond the linear stability analysis to extract any conclusion about 
the stability of the critical point, which can be done by considering specihc models. In 
Section 4, we include some interesting models for which the dynamics of the universe 
is completely independent of the material content. This study is of great relevance 
for understanding the dynamical evolution of more complicated models in the later 
sections. In Section 5, we consider a subfamily of models, which we denote the “term- 
by-term” models for reasons that will become apparent along the article. As we will 
show these models can describe a matter dominated universe followed by a de Sitter 
expansion. In Section 6, we consider what we call the “tripod” models, which are based 
on three potentials. These models can be the simplest examples of models belonging 
to this family which reproduce a consistent cosmological history. Finally, in Section 7, 
we summarize and discuss our results. 
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2 The model 


The most general Horndeski models with a spatially flat de Sitter critical point for 
any material content can be classified in two families [29] . Their main characteristic is 
that the first family has a point-like Lagrangian linear on the held derivatives and the 
second family has a minisuperspace Lagrangian which vanishes on shell in a, that is, 
when it is evaluated at the scale factor of the critical point. The point-like Lagrangian 
for these two families are written as [29] 


L 




I r. 




( 2 . 1 ) 


where 

Teh = —3ilTp2 o,^ , 

is the usual Einstein-Hilbert Lagrangian with H the Hubble parameter, 


.^linear ® E [3VAU,{<P) + ^W,{<P) 




i=0 


( 2 . 2 ) 


(2.3) 


is the Lagrangian of the first family of models, where Ui and Wi being functions of 0, 
and 

3 

^nonlinear = ^ (^0, 0^ H\ (2.4) 

i=0 

describes the second family, with Xi being arbitrary functions. In addition. 


Lm = -a^p{a) = -a^^ps(a), with ps{a) = ps,oa (' 2 . 5 ) 

S 

is the Lagrangian for the minimally coupled matter species labeled by s, and the 
conditions 

3 3 

5 ^ Wi{4>)A‘'^ = ( 2 . 6 ) 

i=0 j=0 


Y,Xi 0 V/" = 0 , 


(2.7) 


have to be satisfied. Conditions (2.6) and (2.7) guarantee that the Lagrangian evalu¬ 
ated at the critical point is linear on 0 which, at the end of the day, is the condition 
needed for the existence of a spatially flat de Sitter critical point for any material 
content (see [29] for technical details). 

In this article, we focus our attention on the first family of models. Therefore, we 
consider that the total Lagrangian is given by 


L Lejj -|- Liinear T L 


( 2 . 8 ) 



(2.9) 


The field equation for these models is given by [29] 




i=0 




H 

7x 






i=0 


Thus, we avoid models for which the left hand side of Eq. (2.9) vanishes, otherwise, 
there is no dynamical evolution of H, as H would also have to vanish. 

The modified Friedmann equation can be expressed using the Hamiltonian density 
which is 

l-i = HeH + "Hlinear + 'Hm = 0, (2-10) 

where the terms are given by 


•Heh = 

3 

•Hlinear = [3(^ - 1) VA t/,(0) + * 0 IT.(0) 

j=0 



( 2 . 11 ) 

( 2 . 12 ) 

(2.13) 


respectively. We have defined for obvious reasons. It must be noted that for the 
model to be able to screen by itself, the Hamiltonian density has to contain at least 
one term with a function of <f>. Thus, we must have at least one IT. 7 ^ 0 with z 7 ^ 0 
for a dynamical approach to a de Sitter vacuum. The modified Friedmann equation 
( 2 . 10 ) can be explicitly written as 

Defining D = 87rGp/{3H‘^) and = 87TGp(j)/{3H‘^), with p^p given by Eq. (2.12), then 
Eq. (2.14) can be expressed as 

D +12^ = 1. (2.15) 

It is also useful to define a dark energy equation of state parameter as 

which should not be confused with the effective equation of state parameter 

u-en = -1 - !§■ (2-17) 


3 Dynamical system 


Defining = Ina, the field equation (2.9) and the modified Friedmann equation (2.14) 
are written as 


H' 


EL H' 






,_/t (1 - ») - 3 ELL" 1) H‘ UiW 


(3.1) 

(3.2) 
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(3.3) 


respectively, subject to the constraint 

i=0 j=0 


where we have dehned Ui = SirGy/AUi/S and Wi = SirGy/AWi/S for notational 
simplicity. In order to close the autonomous system we assume that the different 
matter species can be treated as noninteracting, then considering hi = his we have 






1 + Ws 2 


Elo H' 

yAt/,,^(0)-hJiy00) 


E 




(3.4) 


which simply arises from the conservation equation, already integrated in Eq. 
and equivalent to 


Q. = 




exp [-3(1 + Ws)N ], 


through 




H' 

3(1 + Ws) + 2 — 
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(2.5) 

(3.5) 

(3.6) 


In the hrst place, taking into account condition (2.6) one can note that Eq. (3.1) 
trivially vanishes for H = \/A, as it should by construction. Thus, there is a de Sitter 
critical point with = y/A for any material content. In the second place, from 
Eq. (3.6), one can note that a critical point which exists for any material content (that 
is, any value of Wi) has to be characterized by = 0, where the subscript c means 
evaluation at the critical point. In third place, considering also Eq. (3.2) one concludes 
that we have a de Sitter critical point given by 


3 

{He = y/A, 4> = 4>c, = 0}, with 0c such as ^^(i — 1) Ui{(j)c) 

1=0 


A 

3' 


(3.7) 


The last relation comes from demanding 0' = hi = 0, 0 = 0c and H = y/~A in Eq. (3.2). 
Although the existence of this de Sitter critical point is independent of the particular 
form of the potentials, its position in the phase space depends on them. Thus, if one 
wants the cosmological model to be well-defined close to the critical point, one has to 
require 0c to be real. If for a particular model there is more than one real solution 0c to 
Eq. (3.7), then we have various de Sitter critical points characterized by He = yfA and 
= 0 (but with different 0c’s). On the other hand, the existence of the critical point 
(or points) given by Eq. (3.7) is not precluding the potential existence of additional 
critical points. Thus, in particular models there could be additional critical points of 
the dynamics apart from our de Sitter critical point (or points). 

Assuming now only one material component and considering linear stability the¬ 
ory, the eigenvalues of the Jacobian matrix of the system given by Eqs. (3.1)~(3.4) at 
the critical point (3.7) are 


A = 0, A = —3, A = —3(1 -|- Ws), 


(3.8) 
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where we have simplihed the resulting expression considering condition (2.6). There¬ 
fore, in order to conclude whether this critical point is stable or not, one should go 
beyond linear order. This is not a trivial undertaking without specifying the partic¬ 
ular potentials. Nevertheless, we can consider the 2-dimensional dynamical problems 
dehned in the sections of the phase space given by the critical values to look for some 
hints about the behaviour of the 3-dimensional system. Firstly, if one chooses to re¬ 
strict attention to a vacuum model, i.e. r 2 s,c = 0 , one also obtains that one of the 
two eigenvalues vanishes. Secondly, restricting attention to the 2-dimensional dynam¬ 
ical problem defined in the section of the phase space = a/A one obtains again 
one eigenvalue A = 0. Thirdly, in the section given by 0 = 0c, the 2-dimensional 
Jacobian matrix evaluated at the critical point has non-vanishing eigenvalues, one for 
A = —3(l-t-tCs) and the sign of the other depending on the particular potentials. Thus, 
for the 2 -dimensional dynamical system dehned in this section, if the critical point is 
an attractor for a model with a given value of Wg such that 1 -|- tCs > 0 , it is also 
an attractor for a different value of Wg > — 1 , whereas it is unstable for models with 
< — 1. As we will see studying particular examples, this dependence of the stability 
on the sign of 1 -|- Wg is present in some 3-dimensional models. 

4 Models in which matter does not matter 

Rehecting on the form of the equations of the dynamical system, Eqs. (3.1)-(3.4), one 
notes that the equation for H' is independent of the material content. Nevertheless, 
in order to integrate this equation to obtain H[N) one must substitute the expression 
for which depends on hlg through Eq. (3.2). Thus, even if there is no direct 

dependence on hlg on Eq. (3.1), the dynamics of the universe generically depends on 
the material content which enters into play through the function 4>{N). Nevertheless, 
it can be seen that for some models the dependence of Eq. (3.1) on 0 can disappear 
and, therefore, the evolution of the universe is independent of the matter content for 
those models. 

In this section, we consider some of these models for which any cosmological 
constant or material is completely screened by the held contribution during the whole 
cosmological history even if the self tuning to de Sitter can take place dynamically. 
Firstly, we will show why the WA potential present in the model cannot be Wq if we 
want to approach a de Sitter solution dynamically. Secondly, we will consider models 
with only one Wi and one Uj. Thirdly, as we need to have at least one potential Wj 
with j 7 ^ 0 but it is not necessary to have any Ui, we will include some comments 
about models with two Wi potentials. Finally, we discuss an additional case for which 
the evolution of H is independent of the material content. 

4.1 Wq. exact de Sitter evolution 

If the only hR 7 ^ 0 is Wq, then H ^ = 0 from Eq. (2.12). The model is not able to screen 
dynamically and it has just an exact de Sitter evolution (incompatible with a consistent 
cosmic history) independently of the number of tA’s present. We can understand this 
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better by noting that condition (2.6) and the field equation (2.9) imply 


3 

= (4.1) 

i=0 

and 

ir„ = ^ ^ (4.2) 

i=0 

respectively. These two equations can only be simultaneously satisfied ii H = a/A 
during the whole evolution. 


4.2 One IT), Uj pair 

Assuming that we have only one ITj and one Uj, with i ^ 0, condition (2.6) implies 

IT, = A"? (4.3) 

which gives the same relation between Wi and Thus, Eq. (3.1) can be written as 

Q /n 

h' = — {h^ - 1), with n = j (4.4) 
% 

where h = H/y/X. This differential equation is independent of 0 or hi,,, therefore, one 
can integrate h' and obtain the dynamics of the universe without any knowledge of the 
material content or the particular form of the potentials. It should be emphasized that 
if n = 0, then we are in a case given by Eq. (2.7) and, thus, h' = 0 during the whole 
evolution. Furthermore, one can integrate Eq. (4.4) to obtain 


1 - h-^ 



(4.5) 


which implies that for h > 1 and n > 0 and for h < 1 and n < 0 we have 


h = 


1-Dexp ( —N 


— Ijn 


(4.6) 


whereas for h < 1 and n > 0 and for h > 1 and n < 0 we have 

- —1/n 


h = 


1 + T>exp ( —N 


(4.7) 


where H > 0 is an integration constant. Therefore, for n > 0 one has that the de 
Sitter critical point is a repeller, since for one branch of solutions h evolves from 1 at 
N —)■ —oo to infinity at a finite N, and for the other branch, to h —)■ 0 when N ^ oo. 
On the other hand, for n < 0 we have that the de Sitter critical point is an attractor. 
This attractor is approached from values of h larger than 1 for solutions described by 
Eq. (4.7), and smaller for those given by Eq. (4.6). 






Let us now focus our attention on the models with a de Sitter critical point, that 
is, modes with n < 0, and consider the branch in which h is a decreasing function of 
N. In this case, considering the limit h ^ 1 of the r. h. s. of Eq. (4.4), we obtain 

3 

h' ~ —h for h S> 1, (4.8) 

i 

for any negative n. Thus, far from the attractor the models behave as having 



Wes ~ 1 if i = 1 : 

(stiff matter). 

(4.9) 


Wes ~ 0 if i = 2 : 

(dust). 

(4.10) 



(curvature). 

(4.11) 


Moreover, considering a particular value of n, one can integrate the exact equations for 
h', Eqs. (4.6) and (4.7), to obtain the time evolution of the scale factor. For example, 
for n = — 1 (which is related to some considerations that we will include in Sec. 5), we 
have 


a(t) 


±D + Cexp 


'3VA 


i/3 


(4.12) 


where the ± signs take into account both branches and C > 0. 

On the other hand, even if some models can reproduce a matter dominated epoch 
at early times, it should be emphasized that the material content is not affecting the 
cosmological dynamics at all. Taking into account Eq. (2.15), one can see that the 
contribution of the material content is completely screened by the field through the 
evolution. We have 


= 1 - 3(j - 1) Uj (4.13) 

The held only appears in this equation but not on Eq. (4.4), hence, it has no effect on 
the evolution of the Universe. 

4.3 One fUj, Wj pair 

If we consider a model with only two Wi potentials and vanishing or constant f/j’s, 
condition (2.6) implies that we have 

w = -a"? Wj, (4.14) 

which j ^ i in order to obtain dynamics, as pointed out after Eq. (2.9). In this case 
Eq. (3.1) can be simplihed to 

1 — 

h' = -3h -(4.15) 

Thus, one obtains 

\ {I - h^-^) \= De-^^, (4.16) 
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which approaches 0 when iV —)■ cx), implying that either the first factor vanishes (being 
h zero or infinity large depending on the sign of j) or h —?■ 1. Thus, at least for some 
of these models, the de Sitter critical point is an attractor. 


Considering the regime h ^ 1 in Eq. 

(4.15) one gets 


h' ~ —-h 

J 

for 

* < J, 

(4.17) 

and 




0 1 

1 

d 

for 

i > j. 

(4.18) 


i 


Thus, for the models which may describe a universe with a Hubble parameter decreas¬ 
ing towards the value of the attractor, one could mimic a cosmological phase dominated 
by stiff matter (i = 1), dust (i = 2), or curvature {i = 3) at early times, as explicitly 
shown in Eqs. (4.9)-(4.11). As in the case of Sec. 4.2, the material content is perfectly 
screened by the field during the whole cosmic evolution. 


4.4 Particular cases 

Let us consider models with 


W, = and (4.19) 


with Oi and f3i constant parameters and the factors A have been included for conve¬ 
nience. Thus, condition (2.6) implies this case Eq. (3.1) is again 

independent of 0 and, therefore, on the material content. Indeed we have 


h' 


q ELo K - ^A) 


(4.20) 


which can be integrated specifying the value of the parameters. Nevertheless, it should 
be noted that this case is genuinely different from those previously shown, since now 
h' depends only on h once we have substituted the particular form of the potentials 
and not just considering which are different from zero. Moreover, one expects that 
other particular choices of Wi and Ui can also lead to a differential equation for h 
independent of 0 and, therefore, a cosmological evolution independent of the matter 
content. 


5 Term-by-Term models (TBT) 

5.1 General considerations 

In this section, we investigate a particular family of models whose evolution depends on 
the matter content and, at the same time, can be studied in some generality without 
restricting the particular form of the potentials from the beginning. These models 
are given by assuming that condition (2.6) is satisfied term by term, “Term-by-Term” 
models. That is, we have 

lE,(0) = f/,,^(0)=A-‘/V,,^(0), (5.1) 
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and, thus, we are left with just four independent potentials. Dehning 




(5.2) 


and again h = H/y/X, we can write the field equation and the modihed Friedmann 
equation, expressed as in Eqs. (3.1) and (3.2), as 



(5.3) 


(5.4) 


respectively. As in the general case, these equations together with the conservation 
of matter, i.e., Eq. (3.4), form an autonomous closed system. Therefore, there is a 
critical point for any value of Ws and number of species s at {he = 1, (j)c, ^s,c = 0}, 
where is such that 


v«{4>c) = v-Mc) + 2r3(<#.c) - i 


(5.5) 


As we have already pointed out for the general case, this critical point is not hyperbolic 


and, thus, the linear analysis is not enough to extract any conclusion about its stability. 


Thus, one has to consider particular models to conclude whether this point is stable 
or not. 


In order to consider a particular model which could describe the current evolution 
of our Universe, let us reflect on the form of Eqs. (5.3) and (5.4). As we have already 


emphasized in Sec. 4, for matter to affect the dynamics one needs to have a differential 


equation for h'{N) which depends on 4>{N), with (j){N) given through Eq. (5.4) in terms 


of Q{N). This requires the presence of at least two potentials with a different functional 
form and with non-vanishing derivatives. Moreover, in order to have self-tuning we 
need to have at least one V^ 0 for i ^ 0. On the other hand, one can study how 
the dynamics of the models in the regions where a particular potential dominates in 
Eq. (5.3) taking into account considerations similar to those presented for the models 
in Sec. 4.2. In the present case, however, matter affects the dynamics, since the 
realization of the different regions of domination as well as their duration depend on 
the particular function (j){N) and, therefore, on fl{N). Nevertheless, this approach can 
suggest which models are candidates to describe a consistent background cosmological 
history. When a given potential Vi dominates we would have an expansion similar to 
that given by Eq. (4.4) for n = —1, namely. 


h' ~ —-(h — 1) 


3 


h' 


(5.6) 


h for h S> 1. 


_ 


i 


i 


It must be noted that Vq can never dominate the denominator of Eq. (5.3) since it is 


not present there. If the critical point is an attractor and h{N) decreases towards this 
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attractor, h = 1, we have a dynamics similar to general relativity for early times with 


Weff — 1 if i = 1 : 

(stih matter). 

(5.7) 

2 

1 + WeS ^ WeS ~ 0 if * = 2 I 

1 

(dust). 

(5.8) 

tyeff ^ 2 if * = 3 : 

(curvature). 

(5.9) 


as in the models of Sec. 4.2. Thus, if we want a model compatible with a matter 
dominated phase before accelerated expansion, we should investigate models with the 
potential V 2 and at least one additional potential. The duration of the cosmological 
epochs as well as their order depend, however, of the material content and the particular 
form of the potentials through Eqs. (5.3) and (5.4). 


5.2 Matter dominated universe with a de Sitter attractor 

As we have just discussed, the dominance of the potential V 2 in Eq. (5.3) in regimes 
with a Hubble parameter much larger than the current value, that is h ^ 1 (which is 
equivalent to Tfo ~ \/A), could lead to a cosmological phase of effective matter 

domination, that is 1 + tCefr — 1- Thus, in this subsection we consider a cosmological 
model hlled only with dust with the aim of describing the current and future evolution 
of our Universe. As we are assuming dust we have 

= (5.10) 

where the star denotes an initial condition. Let us consider a model with 

Vo = -10-V and V 2 = (5.11) 

and with the other potentials vanishing. Then Eqs. (5.3) and (5.4) can be written as 

+V) ’ 

_ (1 - + 3 X 10-6 (0 - 

^ ~ 2 X 10-6 

Considering these expressions in Eq. (5.5), one obtains that the model is well-dehned 
close to the critical point, as we have a real value of the held at the critical point, that 
is, (pc = 12.72. One can explicitly check that this critical point is indeed an attractor 
by depicting different sections of the phase space diagram at the critical values. In 
Figs. 1 and 2, we show that this is indeed the case. Moreover, if we consider that the 
universe is hlled with a diherent huid also with 1 + tc > 0, as tc = 1/3, the critical 
point is also an attractor as suggested at the end of Section 3. 

One can numerically integrate the equations of the dynamics, given by Eqs. (5.3) 
and (5.4), to represent the dynamical quantities which caracterize the model. In Fig. 3, 
we show the evolution of the matter parameter Om with respect to 

logi„(l + j) = [(N - N,)-\n(l + z,)] , (5.14) 


(5.12) 

(5.13) 
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Figure 1. Section of the phase space diagram for a model with Fq = ~10 ^md V 2 = 
The diagram depicts a section at the critical value of nm,c = 0. 


w = 0 w = 1/ 3 



10 The diagram on the left and on the right show sections at (j)c of the phase space 

diagram for a model only with dust matter and only with radiation, respectively. 

where 2 ; is the redshift, = z{N^). We also compare the energy densities in logarith¬ 
mic scales, In(Kpni) and \yi{kp^) with k = Sir G/{3A). As can be seen in Fig. 3 the 
solution describes an early matter dominated universe that then evolves to an accel- 
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erated expansion driven by the scalar field. Taking into acconnt the evolntion of the 
effective eqnation of state parameter given by Eq. (2.17), we note that even thongh its 
valne at present is within cnrrent bonnds, the valne of flm is smaller than the desired 
in order to satisfy observational constraints. Thus, even if the model has interesting 
dynamical properties, it does not entirely reproduce the cosmological evolution of our 
Universe. 

One could try to fine tune the parameters of the model or even change the func¬ 
tional form of the potentials in order to obtain a cosmological evolution closer to what 
observations suggest. However, this is not a trivial task because the considered model 
has some successful characteristics due to a long domination of U 2 , which implies a 
long period of the evolution of the universe where the dynamics is simply given by Eq. 
(5.7). This behaviour is also suggested by the fact that a model with 

Fi = -IQ-V and F 2 = (5.15) 

leads to exactly the same dynamics as the one depicted in Fig. 3. We are not suggesting 
that the cosmological behaviour of the models with V 2 and an extra V) to be trivial. 
For example, a model described by 

V 2 = and F 3 = -10“^, (5.16) 

with the same initial conditions, the scalar field quickly dominates the evolution of the 
Universe and does not give way to a matter dominated era. This model is, therefore, 
inviable. However, a change of the initial condition for the field, for example 0* = 10^, 
affects the dynamics in such a way that one recovers the same behaviour as in Fig. 3. 


6 Tripod models 

Although the models studied in the previous section can reproduce a matter dominated 
phase before the accelerating expansion, it is not trivial to reproduce a long enough 
radiation dominated epoch with those potentials. In this section, we investigate the 
simplest models which can, in principle, deliver a viable cosmological history of our 
Universe. To this aim we proceed as in the previous section, studying which kind of 
terms must dominate in order to describe this evolution. As these models are based 
on three nonvanishing potentials we call them the tripod models. 


6.1 General considerations 

Let us write schematically Eq. (3.1) when only one F or IF term dominates in the 
numerator and denominator of each of the terms. Then we have. 


H' ~ 3FA 


iWiHi iWiH^ 


iWi 



( 6 . 1 ) 


Domination of the second term leads to similar conclusions as in Sec. 4.3, that is, we 
can only obtain a dynamical evolution compatible with stiff matter {i = 1 ), matter 
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- log 10 (i+z) 



- log 10 (i+z) 



-log io(i+z) 


Figure 3. The top panel shown the evolution of In(Kpni) and ln(K/?^) for a model with 
Vo = — and V 2 = filled only with dust. The initial conditions used are 

z* = 10^, /i* = 17.84 X 10^, (/>* = 1 (a similar figure can be obtained with = 10^), 
^m* = 1 (solid line) and = 0.99 (dotted line). The middle panel show the evolution of 
flm and Q(j, and the bottom panel tCefr for Qm* = 1- The present values flm, now = 0.315 and 
'a’eff, now = —0.685 are depicted in the middle and bottom panels with dashed lines. 
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domination {i = 2), and curvature domination (i = 3). Thus, we consider that one of 
the terms of our model has to be W 2 - As we want to describe a radiation dominated 
phase similar to that produced in general relativity, we also need a regime during 
which the hrst term also dominates. Considering the hrst term on the r.h.s, we require 
j = i + 1 such that it is proportional to H (as in the dynamic equation for general 
relativity). In addition, the models must contain at least three potentials U and W for 
matter to affect the dynamics (see subsection 4.2). Essentially, if there are only two 
potentials, we can reduce the system to a one potential case by constraint ( 2 . 6 ) and 
consequently the held dependence cancels out in the ratio of the potentials in Eq. (3.1). 
Therefore, we consider that the simplest model must contain the potentials W 2 7 ^ 0, 

7 ^ 0 and one additional potential. As the model with only W 2 and U 2 leads to 
interesting results for matter domination, we focus our attention on the model with 
W 2 7 ^ 0, 7 ^ 0 and Uz,<f> 7 ^ 0. It must be noted that this argument has helped us to 

hnd a particular working model, which does not mean that there are no other tripod 
models that may also be able to describe our Universe by choosing a different third 
potential. 

Taking into account W 2 7 ^ 0, 7 ^ 0 and 7 ^ 0 in Eqs. (3.1) and (3.2) we 

obtain 


h' 


3 

~2 
1 - 


(h-1) 


W 2 ’ 


Q + 3(2h - 1)U2 - 

2hW 


6hf W 2 d(j) — Qhc 
2 


( 6 . 2 ) 

(6.3) 


where U 2 and W 2 , are dehned as before, and we have considered condition ( 2 . 6 ) to 
substitute 


t/3 


A-V2 


IU 2 d 0 -U 2 + C 


(6.4) 


The constant, c, specihes the particular potential t/ 3 . It must be kept in mind that for 
the model to be well dehned close to the critical point (when 0 ' = 0 , hi = tie = 0 and 
h = he = 1), the following relation 


1 + 3t/2 - 6 


W 2 d(j) — 6 c = 0 , 


(6.5) 


must have at least one real solution 0 c, as already emphasized in the general case 
through Eq. (3.7). 

If the de Sitter critical point is an attractor of the model and the Hubble parameter 
decreases towards the critical point, early enough (h S> 1 ) one has 


h' 


3 , ^2,0 
2 W 2 


1 + tCeff ~ 


W 2 ■ 


( 6 . 6 ) 


Therefore, to describe a cosmology compatible with our Universe we need the presence 
of potential terms such that the quotient of these two potentials evolves from 4/3 at 
early A^ to 1 at larger N. As the quotient is a function of 0 we need to consider the 
particular material content to get the evolution in terms of N as required. 
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6.2 Matter and radiation with a de Sitter attractor 

Let us consider a model with 

and W 2 = Ae^'^ + I3e^^. (6.7) 

3 


7/3 is given by Eq. (6.4), and the other potentials vanish. Equations (3.1) and (3.2) 
can then be written as 


h' = 


1 - + 3 (2h - 1) - 6h + e^'^) - Qhc 

2h (Ae^7 + /5e/3<A) ’ 


( 6 . 8 ) 

(6.9) 


respectively, where A > 0 and (3 > 0. The potentials (6.7) have been chosen to 
reproduce an epoch of matter domination and another of radiation domination. This 
can be noted by considering the regime h 3> 1 in Eq. (6.8), which implies 


1 + Weff - 


Xe^'3 


( 6 . 10 ) 


Thus, one has Wes — 0 when the term dominates and tCefr — 1/3 when is 
dominant. The actual occurrence and duration of these phases depend, however, on 
the behaviour of the function 4>{N). If this function is monotonic, then both regimes 
can be attained during the evolution of the universe. Therefore, we will choose the 
value of the parameters A, (3 and c to ensure a positive or negative sign of 0', at least 
while h ^ 1. If we assume that the dark energy contribution is negligible when the 
/3-term dominates, 72 ~ 1, and that this happens when h ^ 1, from Eq. (6.9) we have 


1 — ‘\rp~P^ 1 

0'~- - - ~ —In (3 c-1-y4e^) , (6-11) 

/J fJ 

In the top panel of Fig. 6 we show the actual evolution of 0 for a particular model, 
where A is an integration constant related with the initial condition. As /3 > 0, one 
would have 0' > 0 if > 3c, and 0' < 0 for < 3c. The hrst case implies that 
either there is a minimum value for the held, 0min, such that 0 > 0min, or that c < 0, 
whereas the second case requires both c > 0 and 0 < 0max- Considering now the 
domination of the A—term, 71 ~ 1 and h S> 1 in Eq. (6.9), one obtains 

0'~ ^ (j) c:^^\n{D — 3>cN), (6.12) 

A A 

with D an integration constant. Thus, the sign of 0' is opposite to the sign of c. This is 
compatible with the considerations in the region of /3-domination and does not imply 
any restriction on the value of 0 if 0' > 0. Additionally, the model is well behaved 
close to the critical point if a real 0c exists, that is a solution of 


1 - 6c = 36^“^^ + 2e^‘^0 


(6.13) 
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by demanding 0' = f2 = 0 and h = 1 sX the critical point. As the r.h.s. of this equation 
is always positive for real values of 0 c; we need to consider models with a potential t /3 
given by Eq. (6.4) with c < 1/6. Moreover, we can consider negative values of c and 
restrict our attention to models with 4>'{N) > 0 for h ^ 1. In this case, one needs 
A > 0 in order to have cosmic expansion equivalent to radiation domination before 
matter domination. As the de Sitter critical point must be reached after the matter 
domination era, i.e., when the A term dominates over the 0 term, with 0' = = 0 and 

h = 1 in Eq. (6.9), we have approximately that 

In ■ (614) 

In Figs. 4 and 5, we show the phase space diagram for a particular choice of the 
parameters of the model. We are considering both matter and radiation, therefore, the 
phase space is 4-dimensional and the critical point is characterized by 

{he = 1, = 0, = 0, 0c}, (6.15) 

with 0c given by Eq. (6.13). Thus, in order to understand the attractor nature of the 
critical point, we have depicted different sections of the phase space diagram. 



Figure 4. Section of the phase space diagram for a model with A = 10, (3 = 2 and c = —0.01. 
The diagram shows a section at the critical values flui.c = 0 and flr,c = 0. 

One can numerically integrate the equations of the dynamical system ( 6 . 8 ) and 
(6.9) considering a universe hlled with non-relativistic matter and radiation. Assuming 
that this happens at large enough values of h and imposing the general relativistic 
value ryefr,eq = 1 / 6 , one can obtain the condition for the held at the matter-radiation 
equivalence, which is given by 

A^‘‘4 (6-16) 
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w = 0 


w = 1/ 3 



Figure 5. Sections of the phase space diagram for a model with A = 10, (3 = 2 and c = —0.01. 
The diagram on the left shows a section at ((c = —0.3584 and = 0 and that on the right 
is depicted using ((c = —0.3584 and nm,c = 0. 


In Figs. 6 and 7 we show the dynamical quantities obtained via this integration in 
terms of —log^g(l + z) = [(Ai — Ngq) — In (1 + Zeq)] /InlO. The central panel of Fig. 6 
depicts the comparison of the evolution of the Hubble parameter of our model with 
that of a ACDMR model (cosmological constant, cold dark matter and radiation) in 
logarithmic scales. Note that our model can reproduce the Hubble parameter of the 
mentioned general relativistic model with great accuracy. As the luminosity distance is 
calculated integrating H{z), this result suggests that our model could equally reproduce 
the expected form of this distance. Moreover, the bottom panel of Fig. 6 depicts 
the evolution of the different energy densities in logarithmic scale. The dark energy 
component tracks radiation during radiation and matter domination until it approaches 
the constant value of the critical point. 

The top panel of Fig. 7 depicts the evolution of the parameters. The qualitative 
behaviour of these parameters is consistent with a viable cosmological evolution. There 
is indeed a radiation dominated era followed by a matter dominated epoch and only at 
low redshift the scalar held takes over the evolution leading the Universe to a de Sitter 
stage. Quantitatively we cannot help noticing that the contribution of the scalar held 
at recombination is certainly larger than allowed by current limits which imply at least 
< 0.02 [30]. In the bottom panel of Fig. 7 we show the evolution of the ehective 
equation of state parameter given by Eq. (2.17). From Fig. 7 we can understand 
that even though the ehective equation of state parameter is within observational 
bounds, the abundance of the scalar held at present, however, is too large. Tuning 
the parameters A, (3 and c, does not improve the dynamics in a substancial form and 
we can only conclude that even though this model presents the adequate qualitative 
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Figure 6. Model with A = 10, /? = 2 and c = —0.01, and conditions imposed at matter- 
radiation equivalence, Zeq = 3.3 X 10^, llni,eq = ^^r,eq = 0.4999, heq = 1.6 X 10^, and cj)eq = 
—0.20118. On the top we show the evolution of the scalar field. In the centre we compare 
the normalized Hubble parameter of our model with that of a ACDMR model with the same 
conditions at A^eq using logarithmic scales. At the bottom we show the evolution of the 
energy density of matter, radiation and the dark energy component. 

















features, it does not pass the observational constraints. 




Figure 7. Model with A = 10, /3 = 2 and c = —0.01, and conditions imposed at matter- 
radiation equivalence, = f^r,eq = 0.4999, heq = 1.6 X 10®, and (^eq = —0.20118. We 

show, on the top panel, the evolution of the parameters of the model and, on the bottom 
panel, the effective equation of state parameter. In dashed lines we include the observationally 
measured values of Hm.now = 0.3175 and tCeff.now = —0.6825. 


7 Summary and further comments 

In this article we have considered the cosmology of the family of models linear in the 
held derivatives of the most general minisuperspace Lagrangian for a scalar held lead¬ 
ing to held equations with no higher than second order derivatives and self tuning to 
a spatially hat de Sitter vacuum. We have studied the dynamical system of these cos¬ 
mologies concluding that the spatially hat de Sitter critical point is not hyperbolic and, 
therefore, one has to consider particular models to conclude whether this point is an 
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attractor going beyond the linear analysis. We have shown that, due to the particular 
form of the Lagrangian the material content affects the cosmological dynamics through 
its effect on the scalar held. We have studied some simple models for which, however, 
matter does not matter, since the differential equation for the Hubble parameter is 
independent of both the material content and of the held value. These toy models 
are of special interest for understanding the behaviour of the diherent potentials when 
they have the dominant contribution in more complex models. 

We have then moved to consider particular models. In the hrst place, we consid¬ 
ered what we called “term-by-term” (TBT) models, since they satisfy the constraint 
for the potentials term by term. As we have seen, it is not difficult to describe an 
early matter dominated universe that then evolves to an accelerated expansion driven 
by the scalar held. Nevertheless, at least for the particular choice of the potentials 
that we considered, the material content decays slightly faster than expected. It must 
be noted, however, that by no means the mechanism is so efficient to immediately 
screen any material content, as suggested in Ref. [28] concerning the models presented 
in Ref. [27]. In the second place, we investigated the tripod models, where only three 
potentials are present, as these models can be the simplest cases capable of producing 
a satisfactory cosmological history. These models are able to describe a radiation dom¬ 
inated epoch preceding matter domination and to have a de Sitter attractor, although 
the material content also decays slightly faster than expected in the particular case 
that we have considered in more detail. Therefore, even though we have not found an 
optimal example, the qualitative good features of these models render them promising 
candidates to describe our Universe’s evolution. 

From the purely formal point of view it is worthy to note that the models con¬ 
sidered in this article may avoid the cosmological constant problem at the cost of 
requiring hne tuning of a number of parameters appearing in the specihc models. The 
cosmological constant problem, however, is not solely based on the hne tuning of a 
constant of nature when it is understood as vacuum energy, but in the necessity of 
hne tuning it depending on the phase transitions that may lead to diherent vacuum 
energies. Thus, the problem could be thought to be alleviated since the hne tuning can 
now be undertaken independently of the underlying particle physics theory describing 
the Universe. 
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